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of spin-parities for a resonances in 22 Mg via resonant elastic scattering of 21 Na+p 

J.J. RVQl.Y. Zhang 1 ' 2 - 3 , A. Parikh 4 - 5 , S.W. Xu\ H. Yamaguchi 6 , D. Kahl 6 , S. Kubono 1 ' 6 ' 11 , J. Hu 1 , P. 
Ma 1 , S.Z. Chen 1 ' 3 , Y. Wakabayashi 7 - 11 , B.H. Sun 8 , H.W. Wang 9 , W.D. Tian 9 , R.F. Chen 1 , B. Guo 10 , T. 
Hashimoto 6 , Y. Togano 11 , S. Hayakawa 6 , T. Teranishi 12 , N. Iwasa 13 , T. Yamada 13 , and T. Komatsubara 14 
1 Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China 
2 School of Nuclear Science and Technology, Lanzhou University, Lanzhou 730000, China 
3 University of Chinese Academy of Sciences, Beijing 100049, China 
4 Department de Fisica i Enginyeria Nuclear, EUETIB, 
Universitat Politecnica de Catalunya, Barcelona E-08036, Spain 
5 Institut d'Estudis Espacials de Catalunya, Barcelona E-08034, Spain 
6 Center for Nuclear Study, University of Tokyo, 
RIKEN campus, Wako, Saitama 351-0198, Japan 
7 Advanced Science Research Center, Japan Atomic Energy Agency, Ibaraki 319-1106, Japan 
^School of Physics and Nuclear Energy Engineering, Beihang University, Beijing 100191, China 
9 Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China 

10 China Institute of Atomic Energy, P.O. Box 275(46), Beijing 102413, China 
11 RIKEN (The Institute of Physical and Chemical Research), Wako, Saitama 351-0198, Japan 
12 Department of Physics, Kyushu University, 6-10-1 Hakozaki, Fukuoka 812-8581, Japan 
13 Department of Physics, University of Tohoku, Miyagi 980-8578, Japan and 
14 Department of Physics, University of Tsukuba, Ibaraki 305-8571, Japan 
(Dated: January 21, 2013) 

The 18 Ne(a,p) 21 Na reaction is thought to be one of the key breakout reactions from the hot 
CNO cycles to the rp-process in type I x-ray bursts. To determine this astrophysical reaction 
rate, the resonance parameters of the compound nucleus 22 Mg have been investigated by measuring 
the resonant elastic scattering of 21 Na+p. An 89 MeV 21 Na radioactive ion beam was produced 
at the CNS Radioactive Ion Beam Separator and bombarded a 8.8 mg/cm 2 thick polyethylene 
target. The recoiled protons were measured at scattering angles of # c . m .«175° and 152° by three 
AE-E silicon telescopes. The excitation function was obtained with a thick-target method over 
energies £4 22 Mg)=5.5-9.2 MeV. The resonance parameters have been determined through an R- 
matrix analysis. For the first time, we have experimentally determined the J^ values for ten states 
above the a threshold in 22 Mg. The 18 Ne(a,p) 21 Na reaction rate has been recalculated, and the 
astrophysical impact of our new rate has been investigated through one-zone postprocessing x-ray 
burst calculations. Our new rate significantly affects the peak nuclear energy generation rate and 
the onset temperature of this breakout reaction in these phenomena. 



PACS numbers: 25.60.-t, 23.50.+Z, 26.50.+X, 27.30.+t 

Type I x-ray bursts (XRBs), one of the most fascinat- 
ing astrophysical phenomena, are characterized by sud- 
den dramatic increases in luminosity of roughly 10-100 s 
in duration, with a total energy release of about 10 39 erg 
per burst. These recurrent phenomena (on timescales of 
hours to days) have been the subject of many observa- 
tional, theoretical and experimental work (for reviews see 
e.g., [il-Q)- The characteristics of XRBs have been sur- 
veyed extensively in a number of space-borne x-ray satel- 
lite observatory missions, including RXTE, BeppoSAX, 
Chandra, HETE-2, and XMM/Newton. More than 90 
galactic XRBs have been identified since their initial dis- 
covery in 1976. These observations have provided abun- 
dant data and opened a new era in x-ray astronomy. The 
bursts have been interpreted as being generated by ther- 
monuclear runaway on the surface of a neutron star that 
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accretes H- and He-rich material from a less evolved com- 
panion star in a close binary system 4, 5j . The accreted 
material burns stably through the hot, /3-limited carbon- 
nitrogen-oxygen (HCNO) [g, lZ| cycles, giving rise to the 
persistent flux. Once critical temperatures and densi- 
ties are achieved, breakout from this region can occur 
through, e.g., a- induced reactions on the waiting point 
nuclei 14 0, 15 and 18 Ne. Through the rapid proton 
capture process (rp-process) [cj-[To|. this eventually re- 
sults in a rapid increase in energy generation (ultimately 
leading to the XRB) and nucleosynthesis up to A~100 
mass region [ll], Qjj]. Among the breakout reactions, the 
18 Ne(a,p) 21 Na reaction may provide the principal break- 
out route from the HCNO cycle into the rp-process [3]; 
however, the actual astrophysical conditions under which 
this occurs depend critically on the actual 18 Ne(a,p) 21 Na 
thermonuclear rate. In the temperature region of interest 
for XRBs, this rate has not been sufficiently well deter- 
mined. 

The reaction rate for 18 Ne(cv,p) 21 Na is dominated by 
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the resonance contributions [13j . The uncertainties of 
current rates are mainly caused by the errors in excita- 
tion energies E x (or resonance energies Er) and reso- 
nance strengths ujj. The first theoretical estimate [l3[ 
of this reaction rate was made based on rather lim- 
ited experimental level-structure information in the com- 
pound nucleus 22 Mg above the a threshold at 8.142 
MeV 14]. After that, the levels in 22 Mg have been ex- 
tensively studied, and more than 40 levels were observed 
above the a threshold. Such high level-density suggests 
that a statistical-model approach might provide a reli- 
able estimate of the rate. However, only natural-parity 
states in 22 Mg can be populated by the 18 Ne+a chan- 
nel, and thus the effective level density will be consid- 
erably lower. The a-unbound states in 22 Mg were pre- 
viously studied by many transfer reaction experiments. 
In the 12 C( 16 0, 6 He) 22 MgJp, 25 Mg( 3 He, 6 He) 22 Mg 
and 24 Mg( 4 He, 6 He) 22 Mg [o] experiments, the excitation 
energies in 22 Mg were determined with a typical uncer- 
tainty of ±20-30 keV. Later on, the excitation energies 
were determined precisely with unprecedented resolution 
by a 24 Mg(p,i) 22 Mg [l8| experiment, in which the uncer- 
tainty was about 1-15 keV for most states above the a 
threshold. With these precise excitation energies, the un- 
certainties in 18 Ne(a,p) 21 Na rate can be largely reduced. 

The above indirect studies mainly focused on the exci- 
tation energies determination, and the spin-parity assign- 
ments were not strictly constrained. Some spin-parity 
assignments were made [l3l . HR [H| simply by referring 
to those of mirror states in 22 Ne; such assignments are 
dubious due to the high level-density in this excitation 
energy region. In a later 24 Mg(p,t) 22 Mg [lj| experiment, 
several spin-parity assignments were made via an angular 
distribution measurement. However, the insufficient res- 
olution of their measurements at scattering angles 8 c . m . 
above 20° makes such ,V assignments questionable [l8| . 
In addition, two new spin-parity assignments were given 
tentatively in our previous low-statistics experiment pol ] , 
and such assignments still need to be confirmed by a high- 
statistics experiment. 

A comparison of all available reaction rates shows dis- 
crepancies of up to several orders of magnitude around 
T~l GK [l8|, and therefore it remains unclear whether 
the statistical-model calculations provide a reliable rate 
estimation in a wide temperature region. At present, 
there are still many resonances (above a threshold) with- 
out firm spin-parity assignments, which need to be deter- 
mined experimentally As a consequence, the accuracy of 
the current 18 Ne(a,p) 21 Na reaction rate is mainly limited 
by the lack of experimental spin-parity and spectroscopic 
information of the resonances in 22 Mg above the a thresh- 
old. 

Up to now, only two direct measurements for the 
18 Ne(a,p) 21 Na reaction have been performed and ex- 
tended down to the energies of E c . m =2.0 MeV [2l| and 
E c . m .=l-7 MeV These energies are still too high 

compared with the energy region -E c . m .<1.5 MeV of in- 
terest for HCNO breakout in XRBs. Exciting results [23[ 



have recently become available at the ISAC II facility at 
TRIUMF, where the 18 Ne(a,po) 21 Na cross section was 
determined in the energy region of E c rn , = 1.19-2. 57 MeV 
by measuring the time-reversal reaction 21 Na(p,a) 18 Ne in 
inverse kinematics. The temperature region of interest in 
XRBs is about 0.4-2.0 GK, corresponding to a Gamow 
energy of E c . m .=0.5-2.8 MeV for this reaction, i.e., to an 
excitation region of £^=8.6-11.0 MeV in 22 Mg. There- 
fore, the recent results [23| are still insufficient for a reli- 
able rate at all temperatures encountered within XRBs. 

In this paper, we determine the thermonuclear 
18 Ne(a,p) 21 Na rate through a new measurement of the 
resonant elastic scattering of 21 Na+p. In this mecha- 
nism, 22 Mg is formed via sub-Coulomb barrier fusion of 
21 Na+p as an excited compound nucleus, whose states 
promptly decay back into 21 Na+p. This process inter- 
feres with Coulomb elastic scattering resulting in a char- 
acteristic resonance pattern observed in the excitation 
function [24{ . With this approach, the excitation function 
was obtained simultaneously in a wide energy range of 
5.5-9.2 MeV in 22 Mg with a well-established thick-target 
method [25T4271 ] . without changing the beam energy. For 
the first time, we have experimentally determined the J™ 
values for ten states above the a threshold in 22 Mg. The 
astrophysical impact on the energy generation rate and 
nucleosynthesis in the XRBs is discussed with our new 
rates. 

In the past ten years, many radioactive ion beams 
(RIBs) have been successfully produced and separated 
at CRIB (CNS Radioactive Ion Beam separator) (28LI29I] 
at the Center for Nuclear Study (CNS), University of 
Tokyo, by using an in-flight method. Most of these RIBs 
were subsequently used in resonant elastic scattering ex- 
periments with a thick-target method 0, [3fl - [33l | , which 
proved to be a successful method as being adopted in the 
present study. For this experiment, some experimental 
details were preliminarily described elsewhere [36j . An 
8.2 MeV/nucleon primary beam of 20 Ne 8+ was acceler- 
ated by an AVF cyclotron (K=79) at RIKEN, and bom- 
barded a liquid nitrogen-cooled D 2 gas target (90 K) [34| 
with an average intensity of 65 pnA. The thickness of D2 
gas was about 2.9 mg/cm 2 at 530 Torr pressure. The 
21 Na beam was produced via the 20 Ne(c?,n) 21 Na reaction 
in inverse kinematics. After the Wien filter, a purity of 
70% for the 21 Na beam was achieved on the secondary 
target. 

The arrangement of detectors and targets at the ex- 
perimental focal plane of CRIB was shown in Fig. 3 
in Ref. [36T Two parallel-plate avalanche counters 
(PPACs) [35| measured the timing and position of the 
incoming beam with a position resolution of 1 mm or 
better. Here, the timing signal was used as part of the 
event trigger and for particle identification, and the po- 
sition and incident angle of each beam particle on the 
target were determined by the two PPACs. The beam 
impinged on an 8.8 mg/cm 2 polyethylene (CH 2 )„ tar- 
get, which was thick enough to stop all the beam ions. 
In addition, a 10 mg/cm 2 thick carbon target was used 
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for evaluating the C background contribution. The 21 Na 
beam bombarded the secondary targets at energy about 
89.4 MeV (AE=1.9 MeV in FWHM). The averaged beam 
intensity was about 2xl0 5 pps. The beam particles were 
clearly identified in an event- by- event mode by using the 
position at PPACb and the TOF between PPACb and 
the RF signal from the cyclotron [36| . The recoiled light 
particles were detected with three Micron [37j silicon AE- 
E telescopes centered at angles of #si=0°, +14° and -14° 
with respect to the beam line, respectively. Each AE-E 
telescope subtended an opening angle of about 10° with 
a solid angle of about 27 msr in the laboratory frame. 
In the center-of-mass (cm.) frame for elastic scattering, 
the relevant averaged scattering angles are determined to 
be c . m .wl75 c , 152° and 151°, respectively. The double- 
sided-strip (16 x 16 strips) AE detectors measured the en- 
ergy, position and timing signals of the particles, and the 
pad E detectors measured their residual energies. The 
recoiled particles were identified by the AE-E method, 
together with a TOF-E method in which timing signals 
were given by PPACb and AE detectors [3^ |. The en- 
ergy calibration for the Si detectors was carried out by 
using secondary proton beams produced with CRIB and 
a standard triple-a source. 

The 21 Na+p elastic-scattering excitation functions 
were reconstructed using the procedure described pre- 
viously [13, IHj]. Figure [1] shows the proton elastic- 
scattering spectrum for a scattering angle of c . m .wl75°. 
The cross-section data were corrected for the stopping 
cross sections of ions in the target [13, HH , and the data 
within the dead-layer region (between AE and E detec- 
tors) were removed from the figure. The carbon back- 
ground contribution, as well as any small contamination 
of high-energy particles in the cocktail beam not rejected 
by the veto E detector, were subtracted accordingly us- 
ing the data taken with the C target. The uncertainties 
shown are mainly of statistical origin. The excitation en- 
ergies indicated on Fig. Q] are calculated by a relation of 
E x =En+Q p . Here, the resonance energy Er is deter- 
mined by an i?-matrix analysis (see below), and the pro- 
ton separation energy Q p =5.504 MeV is adopted based 
on the updated masses of 21 Na and 22 Mg [Jj, H|[. The 
corresponding energies adopted in Ref. [181 ] are also indi- 
cated in the figure for comparison, which agree well with 
ours within the uncertainties. 

The 21 Na+p excitation function has been analyzed by 
a multichannel i?-matrix [|(| code MULTI [IT]. An over- 
all i?-matrix fit is also shown in Fig. [TJ A channel radius 
of J4=1.35(l+213) fm [TKU was adopted in the cal- 
culation. The successful reproduction of the well-known 
states [Till at 6.329, 6.587, 6.611, and 6.792 MeV by 
the code (see Fig. [1]) provides confidence in the present 
method. In this paper, we focus on determining the res- 
onance parameters of those states above the a threshold 
in 22 Mg, which eventually determine the 18 Ne(a,p) 21 Na 
reaction rate. The resonance parameters for all observed 
states will be published elsewhere in more detail [42j |. 

In total, ten resonances above the a-threshold were 
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FIG. 1: (Color online) Experimental cm. differential cross 
section for resonant elastic scattering of 21 Na+p at a scat- 
tering angle of # c . m .~175°. It also shows a best overall R- 
matrix fit. The energies adopted in Ref. [l|| are indicated in 
the parentheses for comparison; the location of a threshold is 
also shown. 



observed and analyzed by the i?-matrix code. For the 
first time, we have experimentally confirmed the J* val- 
ues tentatively assigned by Matic et al. [18] for seven 
states at 8.181, 8.519, 8.574, 8.783, 8.932, 9.080 and 9.157 
MeV, and assigned here new values for three states 
at 8.385, 8.657 and 8.743 MeV. As an example, the typ- 
ical i?-matrix fits with possible J 71 ", channel spin s and 
orbital angular momentum £ for the presently observed 
8.578, 8.353, 8.677 and 8.727 MeV states are shown in 
Fig. The presently observed 8.578 MeV state is clos- 
est to the 8.574 MeV in the work of Matic et al. in which 
it was assigned to be a 4 + state based on a shell-model 
calculation. As shown in Fig. [U(a), both 2 + and 4 + can 
fit our data very well. As such, our data support the pre- 
vious 4 + assignment. The observed 8.353 MeV state is 
regarded as the 8.385 MeV state of Matic et al. whose J™ 
was suggested to be 2 + by referring to the mirror state in 
22 Ne. In addition, we assigned it J 7r =(l + -3 + ) in a pre- 
vious low-statistics experiment [20j where 1 + was also 
the most probable assignment. In this work, J 7T =1 + is 
again the best candidate as shown in Fig.^b). Further- 
more, this state was only weakly populated in the previ- 
ous transfer- reaction experiments [la . ITtI [Hj] which pref- 
erentially populated the natural-parity states in 22 Mg. 
This supports our assignment of 1 + unnatural-parity to 
this state. The observed 8.677 MeV state corresponds to 
Matic et al. 8.657 MeV state, which was assigned as a 
J 7r =0 + also based on a shell-model calculation. However, 
such a prediction is questionable because of the high level 
density at such a high excitation energy region. Matic et 
al. regarded this state as the 8.613 MeV state observed 
in Ref. [l5|] where it was assumed to be 3~ by simply 
shifting the energy of mirror 8.741 MeV state in 22 Ne by 
^130 keV. As shown in Fig. [He), a J 7r =2 + is now ex- 
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FIG. 2: (Color online) Sample i?-matrix fitting results for 
some resonances above a-threshold. The (red) thicker lines 
represent the best fits. The relevant channel spin s and orbital 
angular moment £ values are indicated. 



1000 



100 



: — T— 1 — 




— 1 ■ 1 ■ 








Present / Gorres et ai. [1 3; 

• • ■ Present/Chen eta/. [15] 

Present / Matic et a/. [1 8] 

Present / He eta/. [20] 
Present/ HF [43] 






■ \ 




■ is -> 

i '• 
















! 



0.5 



1.0 

r 9 



1.5 



2.0 



FIG. 3: (Color online) Ratios between present reaction rate 
and previous ones as a function of temperature. 



elusive ly assigned to this state. The observed 8.727 MeV 
state is regarded as the Matic et al. 8.743 MeV state, 
which was simply assumed to be the mirror of the 8.976 
MeV, 4 + state in 22 Ne. The present i?-matrix fit strongly 
prefers a 2 + rather than a 4 + as shown in Fig. ^d). It 
is worth mentioning that the data at scattering angle 
around c . m .~152° also support the J 77 assignments dis- 
cussed above [42| . 

The 18 Ne(a,p) 21 Na rate is calculated by a narrow res- 
onance formalism [l8|]. Part of the resonance pa- 
rameters for the reaction rate calculations are summa- 
rized in Table [H The proton partial widths (IV) de- 
duced from our data will be given elsewhere [42j as 
here we have calculated the 18 Ne(a,p) 21 Na rate assuming 
iL>7= " r ° p ~uT a . In this calculation, all resonance ener- 



gies En and strengths W7 are adopted from the work of 
Matic et ai, except for those states with new J 7 ' values 
whose strengths are recalculated as listed in Table [I] Res- 
onance parameters for states between E x =8. 743-13. 010 
MeV have been adopted directly from Table VII of 
Ref. [HI for the present rate calculation. In order to 



TABLE I: Part of resonance parameters utilized for the 
18 Ne(a,p) 21 Na rate calculation. Resonance parameters for 
states between E x —8. 743-13. 010 MeV have been adopted di- 
rectly from Table VII of Ref. [l8l ] for the present calculation. 
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FIG. 4: Nuclear energy generation rates during one-zone XRB 
calculations using the K04 thermodynamic history 
suits using the present rate (black) and Gorres et al. 
(red) are indicated. 



demonstrate the discrepancies between different rates, 
the ratios between the present rate and the previous ones 
are shown in Fig. [3] Comparing to the rate of Matic et 
at, the present rate is much smaller below 0.13 GK, be- 
cause of the 1 + character that we have now assigned to 
the 8.385 MeV state. Consequently, this state does not 
contribute to the rate due to its unnatural-parity prop- 
erty. On the other hand, the present rate is about 2 
times larger around 0.2 GK, because of our new 2+ as- 
signments for the 8.657 and 8.743 MeV states. Compar- 
ing to other available rates [l3|, UU, [20[ , our new rate is 
about a factor of 2-1000 times larger within the temper- 
ature region of interest for XRBs. The curve labeled HF 
is a Hauser-Feshbach statistical-model calculation taken 
from Ref. 43j. It shows that our new rate and Matic et 
al. rate are about 2-5 times larger than the theoretical 



prediction beyond 1.0 GK. 

The impact of our new 18 Ne(a,p) 21 Na rate was exam- 
ined in the framework of one-zone XRB postprocessing 
calculations. Three different XRB thermodynamic his- 
tories were employed, including the K04 (T pea k=1.4 GK) 
and SOI (T poak =1.9 GK) models 0111, as well as that 
of the hottest shell in burst 3 (Xp C ak=l-3 GK) of the hy- 



5 



drodynamic model 1 from Ref. [46|. For each of these 
histories, separate postprocessing calculations were per- 
formed using the previous rates from Gorres et al. [13], 
Chen et al. HH, Matic et al. [Hj|, He et al. and the 
present rate; all other rates in the reaction network [44J 
were left unchanged. For brevity, we only discuss below 
the impact of the Gorres et al. and present rates in the 
K04 model. For both the K04 and SOI models, the rate 
of Matic et al. gave qualitatively similar results to those 
using the present rate; the rates from Refs. [H, [2(| gave 
similar results to those using the rate of Gorres et al. 
Furthermore, no significant differences in either the nu- 
clear energy generation rate or the final nucleosynthesis 
were observed when comparing results using any of the 
above rates with the profile extracted from Ref. [46| . 

As shown in Fig. 01 a striking difference in the nuclear 
energy generation rate at early times is seen when com- 
paring XRB calculations using the present and Gorres 
et al. rates. The nuclear energy generation rate in- 
creases by a factor of 1.8 between 0.3 to 0.35 s using the 
present rate. This also results in a dramatic change in the 
18 Ne(a,p) 21 Na reaction flux at these times. For example, 
at 0.35 s, the 18 Ne(ev,p) 21 Na reaction flux is increased by 
about a factor of 3. Our relatively larger new rate re- 
sults in the depletion of 15 and 18 Ne at early times by 
roughly a factor of 4 relative to abundances calculated 
using Gorres et al. rate. 



In addition, the 18 Ne(a,p) 21 Na rate dominates over the 
/3-decay of 18 Ne at an onset temperature of T 9 «0.47 (as- 
suming a typical XRB density of 10 6 g/cm 3 ) with the 
Matic et al. and present rates. This critical tempera- 
ture is noticeably lower than the breakout temperature 
of TgwO.60 determined using the rates from the work of 
Gorres et al. and Chen et al. It implies that this reaction 
initiates the breakout earlier than previously thought. 

Given the dramatic impact of the 18 Ne(a,p) 21 Na rate 
on the nuclear energy generation rate in these postpro- 
cessing calculations, it is clear that similar tests using 
full hydrodynamic XRB models are urgently needed to 
examine these effects in detail. Such tests are already 
underway. 
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